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ABSTRACT 

We show that a dual-population unification scheme provides a successful paradigm 
with which to describe the evolution and beaming of all bright extragalactic ra¬ 
dio sources. The paradigm consists of two intrinsic radio-source populations, based 
on the two distinct radio-galaxy morphologies of Fanaroff-Riley classes I and II. 
These represent the ‘unbeamed’ or ‘side-on’ parent populations of steep radio spectra; 
the ‘beamed’ source types including flat-spectrum quasars and BLLac objects, arise 
through the random alignment of their radio-axis to our line-of-sight where Doppler- 
beaming of the relativistic radio jets produces highly anisotropic radio emission. 

We develop the model in two stages. In the first stage the source space density 
as a function of cosmic epoch is determined for the two parent populations, and for 
this we use low-frequency source-count and identification data to avoid biases due to 
Doppler-enhanced radio emission. The second stage defines the beaming models for 
each population, using high-frequency survey data and in particular the 5-GHz source 
count in which at high flux densities the flat- and steep-spectrum sources contribute 
in similar measure. We assume that the flat-spectrum objects, quasars and BL Lac 
objects are ‘beamed’ versions of FRI and FRII objects in which the close alignment of 
radio-axis with line-of-sight has changed the radio appearance into a core-dominated 
(flat-spectrum) object. We adopt a simple parameterisation of the beaming, orient 
the parent populations at random with a Monte-Carlo process, and use a minimisaton 
process to determine beaming parameters which yield a best fit to the 5-GHz souce 
count. The best-fit parameters are found to be in good agreement with those measured 
observationally for individual radio sources. In this, the model accurately reproduces 
the change in source-count form with frequency. Indeed the unified-scheme paradigm 
has great predictive power, and we show how the model successfully describes several 
additional and independent data sets. 

Key words: galaxies:evolution, galaxies:jets, quasars:general, radio contin¬ 

uum: galaxies 


1 INTRODUCTION 


The dramatic space density evolution of radio-loud AGNs 
has been shown to mirror the cosmic star formation rate 


(Wall 1998; Boyle & Terlevich 1998; Shaver et al. 1998). 


galaxy population (Metcalfe et al. 1995; 

Roche et al. 1996 

), 

IRAS galaxies (Saunders et al. 199C 

) and the radio-bright 

starbursting galaxies (Condon 1984; 

Rowan-Robinson et al. 


1993). These correlated epoch-dependences provide enticing 
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clues towards understanding AGN physics, in the context 
of triggers to activity such as coalescence of critical mass, 
merging, fuel supply, and the lobe confinement provided by 
an evolving IGM. 


Our concept of the radio-source p opulations devel¬ 
oped slowly. The work o f |Longair (1966 )| an d collaborators 
(Doroshkevich, Longair & Zeldovich 197C) demonstrated 
conclusively the dependence of the radio-AGN evolution 
on radio luminosity: at redsliifts > 1 the most powerful 
radio sources showed space-density enhancements by fac¬ 
tors of > 10 3 over present-day densities, whereas the less- 
powerful sources_showed_Jittleor no enhancement. Subse¬ 
quently Fanaroff fe Riley (1974) demonstrated a morpho- 
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logical dichotomy: double-lobe radio sources with regions 
of highest surface brightness close to the nucleus (Fanaroff- 
Riley type (FRI)) have radio luminosities consistently lower 
than those with highest surface brightness at the extremities 
(FRII). The FRI - FRII division occurs approximately^ at 
a radio power of log 10 (Pisi mhz ) ~ 10 25 W Hz " 1 sr _1 , close 
to the luminosity dividing strongly-evolving radio sources 
from those showing litt le evolution. On the basis of this co¬ 
incidence, Wall (1980) suggested that the FRI population 
showed near-constant space density while the FRII objects 
evolved strongly with cosmic epoch. 

The early surveys were at relatively low frequencies 
and catalogued predominantly steep-spectrum radio galax¬ 
ies of extended double-lobed structures. During the late 
1960s and the 1970s surveys at frequencies > 2 GHz found 
a large number of radio AGN whose spectra did not de¬ 
crease monotonically with frequency - sources of hard ra¬ 
dio spectra, described (unfortunately) as ‘flat spectrum’ ob¬ 
jects. In contrast to the host objects of steep radio spec¬ 
tra, predominantly giant ellipticals, the majority of the flat- 
spectrum objects were identified with QSOs (quasars) or 
BLLac objects and had compa ct radio structu res, unre- 
solved until mappe d with VLBI. Schmidt (1976) and 


Mas¬ 


son & Wall (1977) claimed that despite their apparent high 


radio luminosity of such objects, they showed little evolu¬ 
tion. This result w as disputed and eventually the compre¬ 
hensive analysis of Dunlop & Peacock (1990) showed that 


qualitatively the flat-spectrum population showed a similar 
evolutionary behaviour to the steep-spectrum population: 
luminosity-dependent density evolution was required. How¬ 
ever, at a particular radio luminosity, the space density evo¬ 
lution of the flat-spectrum population is less than that of the 
steep-spectrum population. At the time, physical connection 
between the steep-spectrum and flat-spectrum populations 
were n ot considered in such analyses. 


optical/UV radiation from the accretion disk system domi¬ 
nates stellar emission. 

If we are to understand the physical significance of the 
evolution which originally suggested different populations, it 
is essential to delineate the populations on a physical basis. 
We cannot map the evolution meaningfully until we know 
the physically-distinct populations for which we wish to de¬ 
rive it. 

The primary purpose of this paper is to advance this 
chicken-and-egg situation in the following way. Firstly we 
consider the dual-population unified scheme, based on the 
premise that the two parent populations are FRI and FRII 
radio galaxies, which give rise to BL Lac objects and quasars 
as beamed progeny. Our basic assumption is that all radio 
sources detected above S$ ghz ~ 3 mjy are encompassed by 
this scheme. Using the best available data in terms of source 
counts and identifications at low frequencies, we estimate the 
space densities of the parent populations. We then set up 
the unified models to define critical angles within which the 
line-of-sight turns these lobe-dominated sources into core- 
dominated sources. Using a Monte Carlo orienting process 
we allow free rein to the fitting process to reproduce a high- 
frequency radio source count. This process results in a model 
for the Doppler beaming of the radio jets with physical pa¬ 
rameters which closely mirror those observed in VLBI mon¬ 
itoring surveys. The dual-population unified scheme makes 
numerous quantitative predictions which can be tested, and 
a number are described in this paper. Crucial further obser¬ 
vations are delineated. 

This is not the first attempt to incorporate a uni¬ 
fication paradigm within sp ace density analyses. Before 


the discovery of dusty tori, Scheuer & Readhead (1979) 


proposed a model of superluminal expansion for quasars 
with a view to explaining the ratio of radio-quiet QSOs 
to radio-loud quasars. A second seminal study by Orr & 


; 1JL UL IU11QJ.U.C1 Cll ill auui dllctiyaco. ___ m lauiu-iuuu L juajaio. i ocguhu ocj.ijlj.hcli ouhcij wi i tx. 

It was the discovery of relativistic expansion (Cohen Browne (1982) proposed that the core-dominated quasars 


et al. 1971; Moffat et al. 1972) through repeated VLBI mea¬ 


were aligned versions of steep-spectrum quasars based on 


surements of the brighter ‘flat-spectrum’ sources which laid 
the foundation for unification paradigms. Whilst relativistic 
motion is the primary mechanism giving rise to anistropic 
radio emission in radio-loud AGN, a se cond mechanism, the 
pres ence of a dusty molecular torus (Antonucci & Miller 
1985) was invoked to explain the optical emission character¬ 
istics of some sources. From these two physical mechanisms 
the paradigm emerged in which the random orientation of a 
‘parent’ population on the plane of the sky gave rise to radio 
sources of apparently diver se type. The unified schemes pro¬ 
posed f or radio-loud AGN (Scheuer fe Readhead 1979 ; Drr & 
Browne 1982; Scheuer 1987^ |Barthel 198S ) have both radio 
and optical appearance governed by proximity of the line- 
of-sight to the radio (presumably rotation) axis. In the radio 
regime, close coincidence changes the morphology from an 
extended steep-spectrum source to a source with a domi¬ 
nant flat-spectrum ‘core’ from the Doppler-boosted base of 
the foreground jet. The optical appearance changes from 
galaxy as viewed side-on to quasar or BLLac object when 
orientation permits a view into the torus opening angle and 


t Note that this division is not a simple function of radio power 
and we return to this point in section 3.4.1. 


the twin relativistic-jet model of Blandford & Rees (1974) 


Using the observed redshift distribution of 3C quasars and 
a simple beaming model they calculated the Lorentz factor, 
7 , from the observed core-to-lobe flux ratio distribution, R, 
for a sample of quasars. They found that the flat-spectrum 
quasar fractions at 2.7 GHz and 5 GHz could be reproduced 
by single values of Rt and 7 . These views of unification were 
limited, primarily in that the dusty torus did not form part 
of the modelling. More recently, Urry and Padovani (1995 
and references therein) considered a simple two-population 
unified scheme, where FRII radio galaxies are the parent 
population of all quasars and FRIs are the parent popula¬ 
tion of all BL Lac-type objects. They fitted the observed 
luminosity function for the steep-spectrum, ‘misaligned’ ob¬ 
jects to the flat-spectrum, ‘aligned’ luminosity function, us¬ 
ing a pure luminosity evolution model coupled with simple 
model of the Doppler beaming which required a large range 
of Lorentz factors. 

In contrast, this paper starts with a definition of the 
two parent populations based on low-frequency data alone, 
and derives independent space-density descriptions for each. 
Using high-frequency data it then derives beaming models 
for each population. A n embryonic version of the present 
model was described by Wall & Jackson (1997); the present 
paper adopts the same thesis and procedure, presents the 
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data in detail, describes a refined version of evolution for 


gair (1983), this dichotomy was reflected by a clean division 


the parent population and develops beaming models which 
provides a distribution of apparent jet velocities as well as 
core-to-lobe flux ratios. In addition, cosmological tests of the 
paradigm are presented. 

Throughout this paper we use h= 0.5 and f2o=l-0. 


2 THE DUAL-POPULATION UNIFIED 
SCHEME 

The current unification paradigm for radio-loud AGN is 
based on two ‘parent’ populations, namely (i) the high- 
radio-power FRII galaxies and (ii) the moderate-radio-power 
FRI galaxies. Both populations exhibit anisotropic radia¬ 
tion ar ising from superluminal motion ('Doppler beaming') 


between high-excitation / strong-emission line objects and 
low-excitation / weak-to-no-emission-line objects. This di¬ 
chotomy in the FRII population is reflected in the ‘dual pop¬ 
ulation’ unified scheme which we propose. We suggest that 
only the class ‘A ’ FRII sources appear as lobe-dominated 
broad-line quasars, or core-dominated broad-line quasars at 
small angles of their radio axes to our line of sight. More¬ 
over, it is only for the class ‘A’ sources that the presence of 
an obscuring torus is invoked. 

The role of the low-excitation FRIIs has not been con¬ 
sidered in unified models or space-density analyses to date. 
However the shape of the radio luminosity function for the 
FRII population shows that there must be many more low- 
excitation sources than high-excitation ones, given that ex¬ 


citation status correlates with radio power (Laing et al 


of the l adio jets. In addition, obscuration by a dusty torus 


1991 


Barthel 1994). If these objects have relativistically- 


contributes to the orientation-dependent appearance of the 
high-power FRIIs; when the torus/rotation axis and the line- 
of-sight are approximately aligned, radiation from the accre¬ 
tion disk about the black hole dominates the galaxy emis- 


beamed radio jets then they must show a beamed sub¬ 
population. There is strong evidence, both direct and cir¬ 
cumstantial, that this sub-population is seen as BLLac ob¬ 
jects. Firstly, around some BL Lacs there are extended struc- 


sion. The orientation angle of the radio axis on the plane 
of the skv determines the observed characteristics and the 

tures consistent in scale and nower with the lobes of FRII 

sources fBrowne et al. 1982; Kolleaard et al. 1992; 

Mur- 

hence the classification of the object. 

phy, Browne & Perley 1993; [Dallacasa et al. 1997). Secondly, 

In the case of the FRII radio galaxy, a side-on view re- 

there are too few FRIs for them to be the only class of host 


veals a lobe-dominated steep-spectrum radio source hosted 
by a giant elliptical galaxy, possibly with narrow emission 
lines of high excitation. As the observer’s line of sight comes 
closer to the radio axis, the object appears as a ‘steep- 
spectrum quasar’, the line-of-sight bringing the blue light 
/ broad-line region of the nuclear accretion disk into view to 
dominate the light of the host galaxy. At close coincidence, 
the radio radiation becomes dominated by the Doppler- 
enhanced core or base of the forward jet; features may show 
superluminal velocities. In the case of the FRI radio galax¬ 
ies, the side-on view shows a galaxy of weak-to-no narrow 
emission lines while the on-axis view reveals a BLLac ob¬ 
ject of essentially featureless optical continuum and compact 
radio structure showing superluminal velocities. 

The steep-spectrum side-on objects, FRI and FRII ra¬ 
dio galaxies plus steep-spectrum quasars, dominate the re¬ 
sults from low-frequency surveys, the complete samples and 
source counts. High-frequency radio surveys favour the ra¬ 
dio sou rces of harder spectrum, i. e. beamed, core-dominated 


galaxy for BL Lacs ( Owen, Lcdlow fc Keel 1996 ). Thirdly 


the lack of observed narrow emission lines in the parent 
sources, in the class-B FRII galaxies, correlates with the 
optical spectral properties of BL Lacs; the majority of radio 
galaxies of intrinsic radio luminosities more powerful than 
t he FRI/FRII divide, show passive elliptical-galaxy spectra 


(Rixon, Wall & Benn 1991). 


Thus our dual-population unified scheme, illustrated in 
Figure 1, consists of two parent populations: (a) the high- 
radio-power FRII radio galaxies which are the parents of 
all radio quasars and some BL Lac-type objects, and (b) 
the moderate-radio-power FRI radio galaxies which are the 
parents of BL Lac-type objects. 

Whilst the identified classes of extragalactic radio 


sources are diverse (see e.g. Wall 1994 for a detailed discus¬ 


sources BLLac objects and core-dominated quasars. Rais¬ 
ing survey frequency selects samples and source-counts dom¬ 
inated by objects with their radio axes close to our line-of- 
sight. 


sion), we adopt the working hypothesis that all radio source 
types found in surveys between 0.15 and 10 GHz down to 
Si .4 ghz ~1 mjy belong to one of the two parent pop ulations 
or to one other population, the starburst galaxies (|Condon 


1984; Windhorst et al. 1985). Table hi describes types of ex¬ 


tragalactic radio source commonly discussed in the li tera¬ 
ture; the only not able exclusion is cluster relic sources ( Gio- 


vanini et al. 1992) which are so few in number and so low in 


The hypothesis that lobe-dominated steep-spectrum 
and core-dominated flat-spectrum quasars are increasingly 
aligned versions of FRII radio galaxies is that originally pro- 


posed by Bcheuer (1987); Peacock (1987); Barthel (1989) 


The correspond ing scheme for FRI - BL Lac unification was 
put forward by Browne (1983); Perez-Fournon & Biermann 


(1984). However, there is accumulated evidence that the 


unified scheme is not quite so straightforward. In particu¬ 
lar the FRII sources are not a homogeneous popul ation in 
terms of their optical/UV spectral characteristics: Hine & 


space density as to have minimal impact on our analysis. 

The FRI and FRII sources are the powerful radio 
sources explicity described by the dual-population unified 
scheme of the previous section. 

The starburst galaxies are late-type, in contrast to 
the other populations in the table. They have intrinsically- 
weak radio e mission, and dominate radio samples at low 
flux densities (Windhorst et al. 1985). There is no evidence 


Longair (1979) showed that FRIIs can have strong narrow 
optical/UV emission lines (class ‘A’) or only weak, if any, 


lines, similar to FRI sources (class ‘B’b Laing_et^iL_£l994] 


showed that in the 3CRR sample of Laing, Riley fe Lon- 


of any Doppler beaming of their radio emission; they lie out¬ 
side the dual-population unified scheme. The contribution of 
this population is included in our analysis using previously- 
determined luminosity functions and evolution. 

The GPS and CSS are ‘peaked-spectrum’ sources, 
objects with peaks in their radio spectra between 0.1 and 
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beamed, flat-spectrum 



Figure la. Unified scheme for FRIT radio sources. 



Figure lb. Unified scheme for FRI radio sources. 
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Table 1. Extragalactic radio source classes. 


Population 

UV/optical 
emission-line 

signature 

Radio 

spectrum 

v ~5 GHz 

Class 

Contribution 

to source count 

V ~1 GHz V ~1 GHz 

FRII, class A 1, 

narrow 

steep 

RG 

prominent 

prominent 


broad 

steep 

RG« 

prominent 

prominent 


narrow 

‘flat’ 

quasar 

< 5 % 

prominent 

FRII, class B 11 

weak/none 

steep 

RG 

prominent 

prominent 


none 

‘flat’ 

BL Lac 

< 5 % 

< 10 % 

FRI 

weak/none 

steep 

RG 

prominent 

prominent 


none 

‘flat’ 

BL Lac 

< 5 % 

< 10 % 

Starburst 

narrow 

steep 

galaxy 

< 1 % 

< 1 % 

GPS 

narrow 

‘flat’ 

RG 

< 1 % 

~5% 


broad 

‘flat’ 

quasar 

< 1 % 

~5% 

CSS 

narrow 

steep 

RG 

~15% 

<1% 


broad 

steep 

quasar 

~5% 

<1% 


b Mine & Longair (1979) emission-line class. 
jt alternatively a steep-spectrum quasar. 


3 GHz 

There is accumulating evidence that GigaHertz- 

(iRobertson 1980; 

Peacock & Gull 1981; 

Dunlop & Peacock 

Peaked 

Spectrum and Compact Steep-Spectrum (GPS and 

199C), advantageous as there is no requirement to ‘guess’ 


CSS) sources represent young stages (< 10 4 yr) of FRII 


radi o so urces (Readhead 1995; Fanti, Vigotti & Di Paolo 
199f |). ‘[ Therefore these source types are encompassed in the 
dual-population unified scheme given that (i) the frequency 
of broad-line obects in these classes is in general agreement 
with the simple model of a broad-line region b eing obscured 
by a surrounding torus (Fanti & Fanti 1990) and (ii) su¬ 


perluminal motion has been observed in nearby GPS galax¬ 
ies, with Lorentz factors similar to those inferred for clas¬ 


sical F RIIs (Giovannini et al. 1995 ; Taylor, Vermeulen fc 
Pearsor 1995). Moreover the shortfall of peaked-spectrum 
quasars may due to Doppler-boosting of the forward radio 
jet flattening the characteristic peaked radio spectra with 
the result that GPS and CSS quasars a re often categorized 
as ‘normal’ FRII quasars (Snellen 1997). 


3 RADIO SOURCE SPACE DENSITIES 
WITHIN THE UNIFIED SCHEME 

3.1 Parametric models of radio source evolution 

The first stage in our analysis is to determine the space den¬ 
sity evolution of the two parent populations, allowing each 
to undergo quite separate evolution histories. For each of 
the two populations we determine its local radio luminosity 
function and adopt a simple parametric description of its 
evolution with epoch. Whilst the parametric approach has 
the disadvantage of applying a rigid form of evolutionary 
behaviour, it has the a dvantage of comprising only a l imited 
set of parame t ers (e. q. WallM^eaxson_&_JjongahU_98C L Or£^ 
Morisawa & Takahara 1987; [Rowan- Robinson 
et al. 1993). The alternative approach is a free-form fit 


the form of the evolution prior to determining an acceptable 
fit. However the free-form method yields a large parameter 
set which is difficult to manage in subsequent analyses and 
which does not necessarily lead to further physical under¬ 
standing. A final consideration in our choice of approach is 
that in describing independent evolution histories for the 
FRI and FRII populations, we are limited to datasets which 
are simply too small to use in a free-form analysis. 

There are two distinct types of parameterized evolu¬ 
tion which may be applicable to extragalactic radio sources. 
In the case of luminosity evolution (LE) the radio luminos¬ 
ity of sources change with epoch. If all sources undergo the 
same degree of LE it is described ‘pure luminosity evolu¬ 
tion’ (PLE). Alternatively in density evolution (DE) the co¬ 
moving space density of a population changes with epoch. 
If all sources undergo the same degree of density evolution 
it is described as ‘pure density evolution’ (PDE). A com¬ 
bination of luminosity and density evolution is of course 
possible. Previous analyses which attributed a single evolu¬ 
tion history to all radio sources (i.e. combining the FRI and 
FRII populations into a single ‘steep-spectrum’ population) 
have concluded that neither PLE no r PDE are applicable 
( e.g. Wall, Pearson fc Longair 198Cl| ) and instead success¬ 
ful descriptions of the data have always required a form of 
luminosity-dependent evolution. 

Furthermore there are two popular forms of parametric 
evolution, applicable to both density and luminosity evolu¬ 
tion, which may have some physical significance: (i) power- 
law evolution of the form (l+z) 71 , (1 + a) -1 being the univer¬ 
sal scale factor, and (ii) exponential evolution of the form 
exp t being the look-back time. Whilst power-law 

evolution was successfully applied to models for optically- 
selected QSOs (Boyle, Shanks fe Peterson 198S), exponential 
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evolution has been found to be more applicable for radio- 
selected samples, as it reproduces the strong cos mic evolu- 
tion ob served at relatively mod est redshift ( e.g. Doroshke- 
vich, Longair & Zeldovich 1970). 


3.2 Methodology 


Our approach follows that of Wall, Pearson & Longair (1980) 
by which the local radio luminosity function is determined 
given a complete luminosity distribution and a chosen evolu¬ 
tion function F. The evolution function F(P, z) modifies the 
local radio luminosity function, po, to give the radio lumi¬ 
nosity function at any epoch z, i.e. p(P, z) = po(P)F(P, z). 
This description is entirely general and can describe lumi¬ 
nosity evolution, density evolution, or a mixture. 

The steps in this process to a description of space den¬ 
sity are: 



24 25 26 27 28 29 


For each population in the source count: 

(i) Compile the luminosity distribution N (P) from a com¬ 
plete flux-limited sample of identifications and redshifts. 

(ii) Choose an evolution function F(P,z). 

(iii) Calculate the local radio luminosity function po(P)'- 


l °9j p ,s, mhz ^ W Hz-’sr-’) 

Figure 2. V/Vmax values (*) for the 137 steep-spectrum FRII 
sources in 3CRR. < V/Vmax > for the total sample is shown 
dashed. < V/Vmax > values for bins in A log 10 (Pi 5 i mHz)= 0.5 
are shown as data points (H - ). 


po{P)dP = 


N(P)dP 

/; (So) F(P,z)dV(z) 


(1) 


so that the radio luminosity function p(P, z) = po{P)F(P, z ) 
is then known at all redshifts. 

(iv) Calculate the contribution to the radio source count: 


r z(S) 


N(> S)= dP 


po(P)F(P, z)dV(z) 


( 2 ) 


where z(S) is the redshift at which a source of intrinsic ra¬ 
dio power P contributes at flux density S. A redshift cut¬ 
off z c may be imposed, beyond which F = 0. In practice 
the differential source count for the population is calculated, 
flux-binned for comparison purposes exactly as the observed 
count has been compiled. 


When the contributions from all populations have been 
found: 


highest radio power FRII sources have undergone far more 
evolution than the lower-power ones. Whilst < V/Vmax. > 
= 0.664 (cr=0.025) for the 137 FRII sources, this value rises 
from 0.415 (cr=0.118) at log 10 (Pisi mhz) = 24.25 to 0.807 
(cr=0.052) at log 10 (Pi 5 iMHz) = 28.25. 

(ii) FRI Evolution. The V/Vmax statistic for the 26 FRI 
sources in 3CRR is shown in Figure In contrast to the 
strong evolution indicated for the FRII population, the FRIs 
show little evidence of evolution, with < V/Vmax > = 0.314 
(cr=0.057) for all 26 sources possibly reflecting some negative 
evolution. The highest-power FRI sources in 3CRR have a 
value of < V/Vmax > = 0.507 ((7=0.144) at log 10 (Pisi mhz) 
= 25.75. 


To model the luminosity dependence of the 
parent-source evolution we adopt exponenti al luminosity- 
dependent-densi ty evolution (c/‘model 4b’ of Wall, Pearson 
& Longair 198C) such that the evolution function F(P, z) 
has the form 


(v) Sum the single-population source counts to give the 
total predicted count. The predicted and observed differen¬ 
tial source counts are compared and the evolution functions 
are adjusted, the process iterating from step (ii) until an 
acceptable fit to the observed count is obtained. 


3.3 Determining the evolution functions F(P,z) 


A suitable trial form for the evolution function can be deter¬ 
mined from the <V/V ma x> statistic for a complete sample. 
The best sample av ailable at low frequencies is t he 3CRR 
sample at 178 MHz ( Laing, Riley fe Longair 1985 and more 
recently published data collated by R. Laing, private com¬ 
munication) as it is complete in terms of radio morphologies 
and redshifts. 

(i) FRII Evolution. A plot ot V/Vmax (Figure^) for the 
137 steep-spectrum FRII sources in 3CRR indicates that the 


F{P, z) = exp M(P)r(z) (3) 

where t{z) is the look-back time in units of the Hubble time. 
For Einstein-de-Sitter (12=1) geometry this is given by 

r(z) = (l-(l+z)- 1 ' 5 ). (4) 


Additionally we apply a redshift cutoff to the populations 
mirroring the observed behaviour of the powerful radio 
sources at high redshift (shaver et al. 1996). The evolution 
function is modified such that the evolution peaks at z c /2 
and declines to zero at the cut-off redshift z c \ 


F = F(P, z) for « < z c / 2, 

F = F(P, z c — z) for Zc /2 < z < z c and 
F = 0 for z > z c . 
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Figure 3. V/Vmax values (*) for the 26 steep-spectrum FRI 
sources in 3CRR. < V/Vmax > for the total sample is shown 
dashed. < V/Vmax > values for bins in A log 10 (Pi 5 i mHz)=0.5 
are shown as data points (H - ). 



Figure 4. Luminosity distributions for the 26 FRI (cross- 
hatched) and 137 FRII (hatched) radio sources in 3CRR. The 
smoothed master luminosity distributions are over-plotted and 
were derived with <7=0.35 in log 10 (Pi5i mHz) as discussed in the 
text. 


The evolution rate M is set between 0 and Mmax dependent 
on radio power P : 

M(P) = for Pi<P< P 2 , 

M(P) = 0 for P < Pi, i.e. no evolution of radio sources 
of radio power less than Pi, 

and M(P) = Mmax for P > P 2 , i.e. sources of radio 
power greater than P 2 undergo maximal evolution. 


3.4 Determining the evolution 

3.4.1 Deriving the LRLF 

Luminosity distributions were compiled for the powerful ra¬ 
dio sources from the complete 3CRR sample, comprising 
173 sources with S 178 mhz > 10.9 Jy. In this analysis the 
10 flat-spectrum sources (oitomhz > —0.5, where S oc v a ) 
are excluded so the remainder are steep-spectrum sources, 
‘uncontaminated’ by the effects of Doppler beaming. 

The 178-MHz flux densities were translated to 151 MHz 
using a single spectral index mhz = —0.75, the mean 
value for the 3CRR steep-spectrum sample. Separate lumi¬ 
nosity distributions were compiled for the FRI and FRII 
steep-spectrum parent populations using the morphologi¬ 
cal classifications compiled by Laing et al. (1983) and by 
R Laing (private communication). The total sample com¬ 
prises 26 FRI and 137 FRII sources. The unbinned data 
were smoothed to give a master luminosity distribution for 
each population by convolving each of the unbinned lumi¬ 
nosities with a Gaussian curve of unit area, varying cr until 
the sum of all contributions is smooth whilst preserving the 
‘real’ features of the distribution. The binned and master 
luminosity distributions for the two populations are shown 
in Figure 
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The local radio luminosity functions for the FRI and 
FRII populations are then determined from equation [l], us¬ 
ing the evolution function described by equation |3| However, 
the scarcity of radio sources in the 3CRR sample with radio 
powers below log | 0 T\ 51 mhz = 10 24 W Hz -1 sr _1 yields an 
incomplete LRLF at these radio powers. This incomplete¬ 
ness is resolved by incorporation of the LRLFs determined 
for (i) the ‘local’ E/S0 galaxies ( Condon 1984 ) considering 
them as part of the FRI population, and fii) the starburst 


galaxies, wit h the LRLF as derived by Rowan-Robinson 
et al. (1993) and the evolution function of [Saunders et al 


(1990). These starburst sources were folded into the space 


density analysis with the set evolution function 


Fstar{z) = exp Qt(z) 

with t ( z ), the look-back time, as given in equation [ij. The 
evolution function F 3 tar describes pure luminosity evolution 
as follows: 


Q = 3.1 for 0 < z < 2, 

Fstar(,Z^ — F s tar(,Z — 2) for 2 Z ^ 5. 
and F sta r{z) = 0 for z > 5. 

The total local radio luminosity function is then the sum 
of the LRLFs for each of the three populations (FRI, FRII 
and starburst galaxies). The derived LRLF (Figure |H[) shows 
that the transition between the FRI and FRII populations is 
gradual and is not a simple transition at some radio power. 
This is also illustrated in Figure 4 where there are a number 
of FRIs in the 3CRR sample with P 151 mhz > 10 25 W Hz ' 1 
sr _1 . The luminosity functions were tapered to avoid dis¬ 
continuities, with the result that the space densities for low- 
power FRIIs (log 10 P 151 mhz < 23.5) and high-power FRIs 
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Table 2. Evolution functions used in the source count fit. 



Population 

Evolution 

function 

Fixed 
evolution ? 

Evolution 

parameters 

FRI 

exp LDDE 

No 

Mmax j 2 c> -^lj P*2 

FRII 

exp LDDE 

No 

Mrnaxt %ci ^1) P2 

starburst 

exp PLE 

Yes 

Q=3.1 


Table 3. Fitted evolution parameter values. 

Evolution parameter values 


Popu¬ 

lation 


Chi-square test 
X 2 ■ ^ 

s'-min 


lo 9 ,„ (P, 51 „„ 2 / W Hz -’sr -’) 

Figure 5. Derived LRLF at 151 MHz. The solid line is the total 
LRLF with constituent populations FRI (dotted), FRII (dashed) 
and starburst galaxies (dot-dashed). 


(log 10 P 151 mhz > 27.5) are insignificant. We find good agree¬ 
ment with other LRLFs determined independently. At the 
low-power en d the total LRLF agrees well with those de¬ 
termined by |Punlop fe Peacock (1990); Sadler, Jenkins & 


FRI Mm,ax —0.0, z c = 5.0 

(Pi & P 2 not used given Mmax=0.0) 

FRII Mmax= 10.93, z c =5.62, 

Pi =25.44, P 2 =27.34 


best fit 


30.73 


33 


Kotanyi (1989); Auriemma et al. (1977). This agreement 

suggests that our FRT LRLF is reasonably well defined. 


degrees of freedom. 


The 151-MHz source count used in this comp arison 
(Figure consists of th e count from the 6 C survey (Hales 


This dt fmition is important as low-power FRI sources arc _ Baldwin fe Warner 1988 ) and the count from the 3CR survey 


far mor e numerous locally than the highor-radio-powor FR.T _ from the 3 CRR catalogue at 178 MHz flLaing, Riley fe Lon- 


found i n the 3CRR sample, and at low flux densities their 
contribution is substantial. 

The model FRI and FRII LRLFs from the Milne (fi=0) 
version o f our successful model fit are in agreement with the 
results of Urry & Padovani (1995), expected as both analyses 


find little or no evolution of the FRI population. However, 
there is a major discrepancy between the two FRII LRLFs. 
Urry and Padovani adopted pure luminosity evolution for 
the FRII population, which predicts a high local space den¬ 
sity of low-power FRIIs. According to Urry and Padovani’s 
LRLF there should be similar local space densities of FRIs 
and FRII s. Much smaller local spa ce densities for FRIIs are 
observed; Ledlow fe Owen (1996) find only 6 % of FRIIs in 
a local survey of Abell clusters (z < 0.09). However, over¬ 
estimation of the LRLF for low-power FRIIs does not af¬ 
fect the analysis performed by LTrry and Padovani which 
concentrates on fitting the RLF for the high-power quasar 
population. 


gair 1983 ), transposed to 151 MHz using a single spectral 
index of -0.75. 

The set of parameters given in Table ^ provides a good 
fit to the observed source count. This fit has strong cosmic 
evolution of the FRII population coupled with no evolution 
of the FRI population. The fit to the differential count is 
shown in Figure j|. For the most powerful FRII sources this 
model predicts a space density enhancement at z= 2.8 of > 
10 4 that of the local space density (Figure Q). That the FRI 
population undergoes little or no evolution is a requirement 
of the 151-MHz source count fit: any significant evolution 
produces an excess of faint sources, in disagreement with 
the decline of the differential source count towards lower flux 
densities. The successful fit also indicates a redsliift cut-off 
in the FRII population, with the fit with z c = 5.62 superior 
to that with z c = 00 at the 99.9% level of significance. As 
a result the model reproduces the ‘quasar epoch’, reflected 
as a peak in the FRII space density (Figure^) around z=2- 
3, as seen for powerful flat-spectrum quasars (|Shaver et al' 


1996). 


3-4-2 Fitting the 151-MHz source count 

The parameters in the evolution functions for FRI and FRII 
sources were optimized by performing a x 2 -minimization be¬ 
tween the observed source count at 151 MHz and the total 
count predicted by the differentially-evolving RLF. The min¬ 
imization was achieved using the AMOEBA downhi ll sim¬ 
plex method in multi-dimensions ( |Press et al. 1 992|) . The 
unique set of evolution parameters thus derived for each par¬ 
ent population is shown in Tabic H. 


4 RADIO-SOURCE BEAMING IN THE 
UNIFIED SCHEME 

According to the dual-population unified scheme, flat- 
spectrum quasars and BL Lac-type sources are the steep- 
spectrum FRII and FRI sources aligned so that radio 
axes are close to our line-of-sight. The result is that the 
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Figure 6 . Source count fit at 151 MHz. Data points represent the observed differential source count at 151 MHz constructed as described 
in the text. The model produces contributions from the three populations, the FRII sources (dashed), the FRI sources (dotted) and 
the starburst galaxies (*). The total model count is shown as a solid line. All counts are shown in relative differential form with 
No = 2400(5i5i mhz ) -1 ' 5 sr -1 . 



Figure 7. Space density enhancements as determined from opti¬ 
mized model parameters for a range of Iog 10 (-Pi 5 i MHz) values as 
shown. 

flat-spectrum emission from the approaching radio jetQ is 

t It is the base of the approaching jet which is flat-spectrum. 
Alternative scenarios are that the relativistic emission is from (i) 
the stationary core of the radio source which is optically thick 
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Doppler-boosted or ‘beamed’. This Doppler boosting can re¬ 
sult in the jet emission dominating the extended lobe emis¬ 
sion. 

Having derived the space density and its epoch- 
dependence for the parent populations, the low-frequency 
counts and identification statistics are satisfied. Translating 
the count of these steep-spectrum objects to 5 GHz reveals 
the shortfall at high flux densities in particular which must 
be there: no flat-spectrum sources have yet been included 
and such sources constitute more than half of those found 
in cm-wavelength surveys. To provide this contribution we 
‘beam’ the parent populations, orienting objects randomly 
to mimic the effects of the Doppler-enhanced radiation when 
lines-of-sight and ejection axes come into close coincidence. 
The parameters required to describe the beamed emission 
for any source are the spectral index of the core/jet emis¬ 
sion, the Lorentz factor 7 = (-\/(l + (v/c ) 2 ) -1 for the bulk 
motion in the jets, and the the intrinsic ratio R c of core- 
to-lobe radio luminosity. We adopt a simple model of the 
Doppler beaming for each of the two populations, with the 
beaming parameters varied to produce the best statistical 
fit to the 5-GHz source count. 


at the observing frequency or (ii) synchrotron plasma emanating 
from the central object. All three scenarios in agreement with the 
results here. 











log i o (R obs ) 
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Figure 8. Observed core-to-extended flux ratios, R 0 b s , for high- 
excitation FRII sources in 3CRR, 47 galaxies (x) and 19 quasars 
(*). Median R 0 b s values for the galaxies, i.e. those sources which 
are unbiased by Doppler beaming are shown as (+). From data 
compiled by R Laing. 


4.1 The Doppler beaming parameters 


We chose to model the Doppler beaming of the FRI and 
FRII populations with a range of intrinsic core-to-extended 
flux ratios (distributed normally about a median value) to¬ 
gether with a single Lorentz factor for each parent popula¬ 
tion. There is an additional complication for the FRI pop¬ 
ulation, for which the intrinsic ratio for the is a function of 
r adio power with significant sca tter about a median value 
(ie Ruiter, Parma & Fanti 199(J. In contrast, there is little 
evidence that the ratio correlates with radio power for the 
FRII population. Figure |i| shows the observed ratio for 47 
high-excitation FRII galaxies from the 3CRR sample. The 
median ratio appears uncorrelated to total radio power, with 
a significant scatte r around the median values. The same re¬ 
sult has been seen ( Morganti et al. 1997) in narrow-line radio 
galaxies from the 2-Jy sample ( [Wall fe Peacock 1985 ). 

For the FRI population we adopt the following function 
from t he obse rved ratio for B2 galaxies (de Ruiter, Parma 
& Fant 1990): 


log 10 Rmed — 0.55 logiQ (-P151 MHz ) T 10.78. 

Where Rmed is the median ratio for FRIs of radio power 
A51 MHz- The intrinsic core-to-extended flux ratio R c is then 
determined as a normal distribution about Rmed with a = 
0.45-R m ed • For the FRII population we reflect the observed 
scatter in the intrinsic ratio R c by simply distributing it nor¬ 
mally about a single median value Rmed, with a = 0.45R me d. 
The derived intrinsic ratio, R c , is frequency-independent, 
apportioning the total source flux density between the ex¬ 
tended and core emission. 


4.2 Modelling the beamed products 


The FRI and FRII populations give rise to beamed and 
unbeamed sources with the observed core-to-extended flux 
ratio, Robs, determining whether a particular source is ob¬ 
served as flat- or steep-spectrum. The procedure is illus¬ 
trated in Figure and described in detail below. Throughout 
all flux densities are K-corrected to observed-frame values 
using a = —0.75. 


(i) A single spectral index, = —0.75 is applied 

to the 151-MHz flux density to determine the flux density of 
the steep-spectrum contribution at 5 GHz. This represents 
the emission from the extended lobes of the source. The 
value of —0.75 is adopted as it is the mean spectral index 
of the stee p-spectrum 3CRR sou rces between 178 MHz and 
750 MHz (Laing & Peacock 1980). We extrapolate this index 
to 5 GHz as there is no evidence that the spectrum of the 
steep-spectrum components in 3 CRR steepens significa ntly 
between 178 MHz and 2.7 GHz (Laing & Peacock 198C). 
This steep-spectrum component flux density, S5 ghz, steep, is 


S5 GHz,steep = Sl51 MHz• (5000/151) °' 75 , 
shown as line A to B in Figure |i| 

(ii) The emission from the flat-spectrum core of the 
source is calculated assuming that the source comprises a 
pair of oppositely-directed relativistic jets of bulk plasma 
velocity (3c. The ejection axis of these jets is aligned at some 
random angle 8 (0° < 8 < 90°) to our line-of-sight. The 
observed core-to-extended flux ratio, Robs, determines the 
degree of beaming observed as follows: 


Robs = Rc A 


where A is the sum of the Doppler enhancement from the 
forward- and counter-jets of the source: 

A = S p f + 51. 

with the Doppler factor for the forward-jet, <5/, as 
Sf = (7(1 — (3 cos #)) -1 
and for the counter-jet, S c , as 
S c = (7(1 + /3cos0 )) -1 
where /3 = v/c = (1 — 7 -2 ) 1 / 2 . 

For radio emission comprising continuously-ejected plasma, 
p = 2 — apiat- We adopt a value for the spectral index of 
the core emission, apiat, of 0.0 as this is the mean spectral 
i ndex almjnfr. of th e core com ponents of the B2 sources 
( |de Ruiter, Parma fe Fanti 1 99(1). 

Therefore the observed core-to-extended flux ratio is given 
by 


Robs = -Rc(( 7(1 - /3cos8))~ 2 + (7(1 + /3 cos 6 »)) —2 ). 

The flux density of the beamed contribution, S5 ghz, flat, 
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is represented by line B to C in Figure |i)] and is simply 
related to the steep-spectrum flux density by Ss ghz = 
Roba-Sb GHz , steep • 


(iii) A source is counted as ‘flat-spectrum’ for values of 
Robs large enough such that a\ ^§Hz — — 0.5. This deter¬ 
mines the value of Rmin, the lowest Robs value at which a 
source is counted as flat-spectrum, as follows: 


The spectral index between 2.7 GHz and 5 GHz is given by 


5 GHz 
»2.7GHz 


log 


!0 S 2 . 


=>5 GHz ,total 
7 GHz ,total 


logi 


For afi a t= 0.0 the flat-spectrum components have the same 
flux density at 2.7 GHz and 5 GHz, i.e. S2.7 Gnzjiat = 
$5 GHz,fiat and for asteep=~ 0.75 the steep-spectrum com¬ 
ponents at 2.7 GHz and 5 GHz are related as S 2.7 ghz, steep 
= 1.6 S's ghz, steep- Thus we determine R m in as having the 
value of 0 . 66 . 


The largest angle between the jet and the line of sight 
for the source to appear flat-spectrum occurs at the the crit¬ 
ical angle 9 C when R a b s = Rmin ■ Thus 9 C = 5 m in = 5(9 C ) so 
that 6 C is determinable from 


cos( 0 c ) 


1 1 / Rc ,p 1 

/3 /3 7 Rmin 


( 5 ) 


assuming that the beaming due to the counter-jet is negli¬ 
gible. In fact our successful model finds that 9 C < 10° and 
the assumption is justified. 

Thus the only unknown parameters in the 5-GHz source 
count fit are values for 7 and R m .ed for each of the FRI and 
FRII populations. 


4.3 Fitting the 5-GHz source count 


The 5-GHz source count is constructed from survey data as 
detailed in Table [i] and is well defined across a wide flux- 
density range. At 5 GHz, the beamed products of the FRI 
and FRII sources make a highly significant contribution to 
the source count; for flux den sities above 0.5 Jy, some 55% 
of sources are ‘flat-spectrum’ (Pauliny-Toth et al. 1978). 

To define the populations in accordance with our ear¬ 
lier discussion of the FRII beamed products, we split the 
FRII parent sources into high- and low-excitation types us¬ 
ing asimglejineai^junction ofjog la (A 5 ]_MHz) as suggested 
by Laing et al. (1994) and Barthel (1994)| . The fraction of 
low-excitation FRIIs was set at 50% at log 10 P 151 mhz = 
25.0, declining to zero at log 10 P 151 mhz = 27.0). 

The source count fit was carried out as for the 151- 
MHz count, this time incorporating the beamed products of 
the FRI and FRII sources by randomly aligning the sources 
with respect to our line-of-sight. The source alignment an¬ 
gle, 9, for each of the FRI and FRII contributions was gener¬ 
ated using NAG pseudo-random number routines. The best- 
fit beaming parameters were determined again using the 
AMOEBA downhill simplex method, minimizing fi 2 eval¬ 
uated between the observed and model source counts at 
5 GHz. 

A good fit to the 5-GHz source count was found for the 
parameter values given in Table 

The observed source count and the count calculated 
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Table 6. Fitted beaming parameter values. 

Chi-square test 


Population 

parameter values % 2 

id 

FRII 

7=8.5, R med =0.01 

9c(Rmed) = 7?1 


FRI 

7=15.0 

p „ p—0.55 

VLmed r 151 MHz 



best fit 32.98 

25 


1 degrees of freedom. 


for the optimum fit are shown in Figure [!(]. Here the sub¬ 
populations contributing to the count are split into their 
beamed and unbeamed products as detailed in Table The 
contribution from the beamed sources peaks at high flux 
densities (0.1 < Ssghz < 10 Jy), broadening the ‘evolu¬ 
tion bulge’ of the count. Thus the unification scheme pro¬ 
vides a straightforward explanation of why the source-count 
‘plateau’ widens with survey frequency. It is the Doppler- 
beaming of the FRII population which has the primary im¬ 
pact in this respect. The combination of zero evolution with 
milder Doppler beaming for the FRI population results in 
a monotonically-decreasing count of beamed and unbeamed 
FRIs throughout the flux density range as shown. 


5 TESTING THE MODEL 

The substantial predictive power of the model provides di¬ 
rect tests of the dual-population unified scheme. 


5.1 The beaming models 

5.1.1 Jet speeds; values of the Lorentz factor 7 


The single most successful aspect of the model is that val¬ 
ues of 7 determined from the optimization process agree 
with those determined from estimates made from (i) VLBI 
observations of superluminal sources and (ii) synchrotron 
self-Compton models for the observed and predicted X-ray 
flux. 


VLBI studies of superluminal sources measure the ap¬ 
parent velocities of recognisable features in radio sources. 
Results from a large (> 100) sample of superluminal sources 
show that, assuming that the bulk and pattern speeds are 
the same, 7 m in ~ 4 are commonly found although there is 
a wide spread of values from 2 < 7 < 20 ( Vermeulen 199E| ). 

The synchrotron self-Compton model estimates the 
Doppler beaming from the observed and predicted X-ray 
flux densities of radio sources. For a sample of ~100 sources 
a lower limit for the bulk Lorentz factor of ~10 is found 
with a viewing angle ~ 8 ° for all superluminal sources. This 
analysis also indicates t hat the core-to-extend ed flux ratio 
has a significant scatter (Ghisellini et al. 1993). 

In addition to the above, 7 ~ 10 successfully accounts 
for the non-detection of a counter-jet in highly-alig ned ob¬ 
jects such as 3C273 (Davis, Unwin fc Muxlow 1991). 
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151 MHz 5 GHz log Frequency 


Figure 9. Modelling the flat- and steep-spectrum components of the FRI and FRII sources at 5 GHz. The contribution from the steep- 
spectrum component at 151 MHz (A) is transposed to 5 GHz (B) using = —0.75 (solid line). If the radio axis lies close to the 

line-of-sight (0 < 0 C ) the core emission (dotted) dominates and the source appears ‘flat-spectrum’; the total flux density of the source is 
*^5 GHz,steep T S 5GHz, flat- 


Table 4. 5-GHz source count data. 

Flux density Survey Survey 

range S 5 GHz instrument(s) name Reference 


10 Jy - 100 Jy 
1.5 Jy - 10 Jy 
0.5 Jy - 1.3 Jy 
67 mjy - 0.5 Jy 
10 mjy - 55 mjy 
1.5 mjy - 8.5 mjy 
1.36 mjy - 6.00 mjy 
16 /rjy - 60 /rJy 


NRAO + MPIfR 

strong sources 

Kuhr, Paulinv-Toth &; Nauber (1981) 


Green Bank 91m 

87 GB 

Gregory &; Condon (1991J 


NRAO + MPIfR. 

S4 

Paulinv-Toth et al. (19781 


NRAO 300ft 

6cm survey 

Davis (19711 



MPIfR. 

deep selected regions 

Panlinv-Totl 

i. Stenne fc Witzel (1980) 


VLA 

E/S0 galaxy 

Wrobel fc Krause (19901 


VLA 

Lynx-2 area 

Donnellv. Partridge V Windhorst (1987) 

VLA 

DEEPS2 

Fomalont et al. (1991) 


5.1.2 The observed core-flux : extended-flux ratio 


The observed core-to-extended flux ratios, Robs, span a 
wide range of values. The most lobe-dominated sources, 
those unaffected by Doppler beaming, have ratios a s low 
5 (e.ff. OD-159 from the 2-Jy sample. 


aS Robs 


10 " 


ganti, Killeen & Tadhunter 1993). In contrast the most 


Mor- 


core-dominated sources are heavily Dop pler-beamed and can 
have R„bs ~ 10 3 ( e.g. PKS B0400+258, Murphy, Browne & 


Perley 1993[ ). 

Our assumption of a normal distribution for R c about 
Rmed results in a wide range of R 0 b s for both FRI and 
FRII beamed products (Figures [T^j and [Ti]), and provides 
a smooth transition between the parent objects and their 
beamed sources. The Robs values for BL Lac-type sources are 
predicted to extend to higher values due to the large range 
of intrinsic R c values. The limited data in complete samples 
appropriate for comparison are in approximate agreement. 


5.2 Multi-frequency source counts 

Figure [k^ shows the observed differential source counts at a 
wide range of radio frequencies together with the predictions 
of counts from the model. The model successfully reproduces 
these source counts, qualitatively at least. There are quati- 
tative differences, most notably the excess in the prediction 
at low flux densities for 1.4 GHz. The most proabale ex¬ 
planation is that our simple model of the local luminosity 
function overestimates the space density for powers below 
log 10 (Pl 5 l MHz) ~ 24.0. 


5.3 Core-dominated and broad-line fractions 

Figure |l4| shows the population mix as predicted by the 
model as a function of 5-GHz flux density. The striking 
feature is the change in dominant population moving from 
high flux densities (high-excitation FRII radio galaxies and 
quasars) to low flux densities (FRI radio galaxies, BL Lac- 
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Table 5. Radio source types at 5 GHz. 


Radio 

Source population/evolution/ Beamed spectrum 

Re$ beamed & unbeamed products at 5 GHz ? a 2 7 GHz 


FRII sources 

tapered exp LDDE, F(P,z) = exp M(P)r(z) 
evolution parameters as determined in section 3 



i 

High-excitation radio galaxies & quasars (class A) 

no 

steep 

2 

Quasars 

yes 

‘flat’ 

3 

Low-excitation radio galaxies (class B) 

no 

steep 

4 

BL-Lac type objects 

yes 

‘flat’ 

FRI sources 

non-evolving as determined in section 3 



5 

FRI radio galaxies (class B) 

no 

steep 

6 

BL-Lac type objects 

yes 

‘flat’ 

low-power 

exp PLE, F(z) = exp Qt{z) 

PLE as described in section 3 



7 

starburst galaxies 

no 

steep 


jj These numbers are used to reference the source types in 


Figures ^ to \v\ . 



Figure 10. Source count fit at 5 GHz. Data points represent the differential source count at 5 GHz from data given in Table |f|. The 
model produces contributions from seven source types as described in Table All counts are shown in relative differential form with 
Nq = 60(S 5 GHz ) -1 ' 5 sr" 1 . 
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z 

< 

\ 

z 

< 



% 10 (flux density / Jy) 


Figure 13. Observed source counts at six frequencies with data poi nts from surveys as given in W all (1994) and the addition of the 
VLA FIRST survey source-count points 1 mjy < Si. 4 GHz A 1 Jy (Becker, White & Helfand 1995). The predicted model counts are 
shown dashed. All counts are in relative differential form A N/ A Nq , where AN is the number of sources per sterad with flux density S v 
between S 2 and .S' 1 and No = K v S’ ,7 1 ' 5 , the number of sources expected in a uniformly-filled Euclidean universe. The values of K v are 
2400, 2730, 3618, 4247, 5677 and 3738 respectively for the six frequencies shown. The horizontal range bars show the flux-density bin 
width S 2 to Si, and the error bars represent y/N uncertainties. Polygons show the count estimates from P(D) (background-deflection) 
analyses. 


type sources and - eventually - starburst galaxies). The dom¬ 
ination of the high-excitation FRII sources at high flux den¬ 
sities produces a pronounced peak in the quasar population 
between 0.5 and 10 Jy; the decline in quasar fraction towards 
lower flux densities is in clo se agreement with obse rvation 
(Figure as discussed by |Wall fe Jackson (1997)|. 


The change in fraction of broad-line sources as a func¬ 
tion of flux density agrees with the best available deter mina- 
tions at 408 MHz (Figure 5 of Wall & Jackson (1997). Such 


a change is evidence for the unified s cheme rather t han evi¬ 
dence against it as was suggested by Singal (1996) . Despite 
its relative simplicity, the dual-population scheme shows 
that tests of unification dealing with broad-brush propor¬ 
tions of broad-line objects or flat-spectrum objects are far 


too simplistic; the mix of sub-populations selected in each 
sample must be considered. There must be similar reserva- 
tions over elementa ry tests based on size discrimination ( e.g. 
Kapahi et al. 1996). 


5.4 The redshift distribution of fiat-spectrum 
quasars 

A single sub-population redshift distribution, N(z), provides 
a potent test of th e model. The flat-spectrum sample of 
Shaver et al. (1996) comprises 444 sources, all of which are 
identified. Of these, 358 quasar redshifts were available to 
us. The sample is moderately complete to S i 2 . 7 GHz= 0.25 Jy 
although a substantial proportion of the sample is limited to 
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Figure 11. Model R 0 bs distribution for the flat-spectrum Figure 14. The predicted population mix at 5 GHz, with popu- 

BL-Lac type sources (cross-hatched), plus steep-spectrum FRI lations as given in Table 

and low-excitation FRII sources (hatched), over 1 sterad with 
§5 GHz >0.1 Jy. 



q W//X//A///Y//X//A///Y//X//A///r//A2QQWjQQ&yjQQQwQysiTTTrn\ 

-3 -2-10 1 

log w (Robs) 


Figure 12. Model R a i, 3 distribution for (i) flat-spectrum quasars 
(hatched), (ii) steep-spectrum quasars and broad-line radio galax¬ 
ies, i.e. high-excitation FRIIs with radio axes within 50° of the 
line-of-sight (hatched), and (iii) narrow-line radio galaxies, high- 
excitation FRIIs with radio axes lying > 50° from the line-of-sight 
(solid fill), over 1 sterad and with SsgHz >0.1 Jy. 


5 *2.7 GHz— 0.6 Jy. All sources were selected as flat-spectrum, 
with al.?GHz > “0.4 (5 ex u a ). The redshift distribution for 
the observed sample of 358 quasars is shown in Figure [l6| 
along with the model prediction for an equivalent complete 
flux-limited sample (52.7 GHz > 0.5 Jy). The model success¬ 
fully replicates the observed N(z) for the largest sample of 
spectroscopically-identified quasars available. It should be 
emphasized that these data are completely independent of 
the (low-frequency) data used to construct the model. 
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Figure 15. Core-dominated fractions at 2.7 GHz. The dashed 
line is the model prediction of core-dominated sources (BL Lacs 
plus quasars). The data points are derived from two sam¬ 
ples discussed in detail by Wall & Jackson (1997), section 3.1: 
PKSCAT90 (+) with S 2 . 7 G.Hz > 0.25 Jy and PSR (o) with 0.10 
< S2.7GHZ < 0.25 Jy. The error in fc is y fiN q ) / (bin total). 


5.5 Radio surveys to mJy flux-density levels 

The new, large-area radio surveys (i.e FIRST, NVSS, 
WENSS & SUMSS) reach mJy flux-density levels and will 
yield samples forming the basis of future direct tests of the 
dual-population model. 

Figure [l7| shows the predicted integral population mix 
as a function of 1.4-GHz flux density. The dominant pop¬ 
ulations in the range 1 < 5i.4GHz < 100 mJy differ sub¬ 
stantially from those at higher flux density. A direct conse¬ 
quence of the dual-population model with its strong differ- 
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Figure 16. The redshift distribution for 358 flat-spectrum 
quasars from Shaver et al. (1996), shown with the model pre¬ 
diction for 358 quasars selected with ?GHz — — 0-4 and 
32.7 GHz > 0-5 Jy. 
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Figure 17. The predicted (integral) population mix at 1.4 GHz 
with the populations as described in Table FI 


ential evolution is the dramatic decline in prominence of the 
high-excitation FRIIs (and hence broad-line objects) toward 
1 mjy. 

This change in population with flux density likewise has 
major influence on the redshift distribution of sources, as 
Table ^ illustrates. At 0.1 < Si. 4 ghz < 100 Jy the popula¬ 
tion is predominantly FRII sources at high redshift, while 
a very much more local distribution of sources appears at 1 
< Si. 4 ghz < 100 mjy as the low-power FRIs and eventually 
the starbur st g alaxies rise to prominence. 

Figure hSl shows a comparison between N(z) distribu¬ 
tions for a flux density range 1 < Si. 4 ghz < 100 mjy. The 
striking differences between our model prediction and the 



Redshift, 2 

Figure 18. The predicted N(z) (per O.Olz per sterad) for sources 
with 1 < Si .4 GHz < 100 mjy. Models are from Dunlop & Peacock 
(1990) (average of their models 1-7), shown as thin lines, and this 
paper, shown as thick lines. Legend: ... flat-spectrum sources, - - 
- steep-spectrum sources and solid line total of all sources. 


average N(z) model of Dunlop & Peacock (1990) 
the sharp spike at z 


(i) 


0.1 due to the starburst galaxies 
in our model, a population not included in the Dunlop & 
Peacock analysis, (ii) the lower median redshift distribution 
of sources in the prediction of the current model, due to 
the dominance of the unevolving FRI sources within this 
flux density range, and (in) the greater prominance of flat- 
spectrum sources in our model, the beamed products of the 
low-power FRI sources. 

Although the model predicts unequivocal N(z) and 
population mix, e.g. for flux densities in the range 1 < 
Siaghz < 100 mjy, it has been developed from consider¬ 
ation of total source counts at high and low frequencies, but 
population mix only at low frequencies and high flux densi¬ 
ties, Sm mhz > 10.9 Jy. Complete samples (with redshifts) 
at the mjy level offer substantial scope for improvement or 
refinement of the par ameters. For example, the sample at 
S 2.7 GHz = 0.3 mjy of Gruppioni et al. (1997) contains more 


broad-line objects than predicted. If confirmed, this suggests 
that we have overestimated the the FRI LRLF at low pow¬ 
ers, a possibility already suggested in §5.2. Tapering off the 
FRI LRLF more rapidly at its higher radio powers would re¬ 
quire a compensatory broadening of the FRII LRLF, result¬ 
ing in an increased proportion of FRII objects at sub-mjy 
levels. 
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Table 7. Population mix at 1.4 GHz. 



FRII high-excitation 

FRII low-excitation 

— 

FRI — 

Starburst 

Flux density range 

BL & NLRGs 

Quasars 

RGs 

BL Lac 

RGs 

BL Lac 

galaxies 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Si .4 GHz > 100 mjy 

59% 

15% 

10% 

4% 

7% 

4% 

1% 

1 < Si. 4 GHz < 100 mjy 

7% 

1% 

4% 

- 

38% 

30% 

20% 


6 CONCLUSIONS 


We have described a paradigm in which evolution of ex- 
tragalactic radio sources is delineated for populations which 
are physically related through orientation-dependent effects. 
The dual-population scheme is based on the premise that the 
two parent populations are radio galaxies of FRI and FRII 
morphologies, with distinctly different evolutions which we 
have determined from analysis of counts and identification 
data at low frequencies. The fundamental assumption is that 
all FRII and FRI radio galaxies have radio lobes fed with rel¬ 
ativistic beams. The beaming models by which these parent 
objects appear as the quasars or BL Lac objects when lines- 
of-sight coincide closely with radio axes were then developed 
by considering the high-frequency counts and identifications, 
dominated by these beamed products. 

The model provides agreement with source-counts and 
with identification and redshift data at high flux densities 
and at both high and low frequencies. It provides a natural 
explanation of the change in shape of source counts with fre¬ 
quency. Furthermore, such a unified scheme provides a nat¬ 
ural explanation of why steep-spectrum sources show more 
pronounced cosmological evolution than do bea med (flat- 
spectru m) sou rces of the same (appare nt) power ( Duulo] 3 _& 


(1979)[ 


Peacocl: 1990), a result anticipated by Scheuer & Readhead 


One of the clearest successes of the model is that the 
beam speeds derived from statistics of counts and identifi¬ 
cations are in reasonable agreement with those determined 
for individual sources from VLBI measurements. Further it 
predicts a core-dominated fraction of sources as a function 
of flux density in excellent agreement with observation. 

Four further remarks are important with respect to the 
current model: 


(i) Although we chose LDDE to describe the evolution, 
it can be shown that due to the shape of the RLF some 
models of LDDE are exactly equivalent to PLE. Obviously 
a simple solution of PLE would be very interesting with 
regard to physical models for all AGN evolution (radio-loud 
and radio-quiet). 

(ii) The model fit of LDDE has the lower-radio power 
FRIIs undergoing little or no evolution. By virtue of the 
correlation between line-strength and radio luminosity, these 
are inevitably of weak emission-line type. This suggests that 
the weak-lined FRIIs underwent a similar evolution history 
as FRIs; it is only the powerful strong-emission line FRIIs 
which evolved strongly with cosmic epoch. It might then 
be possible to define the two populations in terms of op¬ 
tical emission-line strength rather than radio morphology, 
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the two emission-types perhaps representing different accre¬ 
tion types. Such a view has important ramifications for the 
process of physical evolution; and it may be that develop¬ 
ment of such a variant of the model could result in a single 
source function describing the entirety of the radio-source 
population. 

(iii) Although the simple beaming models adopted pro¬ 
vide a good description of the count data, there are other 
combinations of the Doppler beaming parameters which 
work. For example a very wide distribution of Lorentz fac¬ 
tors also fits the GHz-frequency source counts and this vari¬ 
ant was adopted by Urry and Padovani (1995). Ultimately 
our selection of parameters was made on the basis that the 
intrinsic core-to-extended flux ratio is an observed quantity 
and that the model distribution of Robs values produced has 
the observed shape. In addition, the single 7 value fitted in 
our model can be interpreted as modelling the mean value 
for the fastest material in the jets. The critical angles in¬ 
ferred from the successful models are ~ 8 ° and ~5° for the 
FRI and FRII populations respectively. 

(iv) It must be borne in mind that the current model has 
been developed fundamentally with only three data sets: a 
source count at 151 MHz, a source count at 5 GHz, and 
a complete set of identifications and redshifts at 178 MHz. 
It is clear that the model can be improved. For example it 
predicts a low flux-density end of the 1.4-GHz count which 
is too high, and it predicts perhaps too few broad-line ob¬ 
jects at low flux densities. Both deficiencies could be recti¬ 
fied with a local radio luminosity function for which the FRI 
component rolls off at somewhat lower powers than the cur¬ 
rent one, compensated by extending the FRII component to 
these lower powers. 

However in view of the limited data from which the 
scheme was developed, its general success in describing much 
of the remaining data provides indication that the concept 
is substantially correct. 


REFERENCES 

Antonucci R., Miller J., 1985, ApJ, 297, 621 
Auriemma C., Perola G. C., Ekers R., Fanti R., Lari C., 

Jaffe W. J., Ulrich M.-H., 1977, AA, 57, 41 
Barthel P. D., 1989, ApJ, 336, 606 
Barthel P. D., 1994, in Bicknell G. V., Dopita M. A., 

Quinn P. J., eds, The First Stromlo Symposium: The 
Physics of Active Galaxies. ASP Conf Ser 54, p. 175 
Becker R. H., White R. L., Helfand D. J., 1995, Ap J, 450, 559 
Blandford R., Rees M., 1974, MNRAS, 169, 395 
Boyle B. J., Terlevich R. J., 1998, MNRAS, 293, L49 














18 C.A. Jackson & J.V. Wall 


Boyle B. J., Shanks T., Peterson B. A., 1988, MNRAS, 235, 935 
Browne I. W. A., Orr M. J. L., Davis R. J., Foley A., 1982, 
MNRAS, 198, 673 

Browne I. W. A., 1983, MNRAS, 204, L23 

Cohen M. H., Cannon W., Purcell G. H., Shaffer D. B., 

Broderick J. J., Kellermann K. I., Jauncey D. L., 1971, ApJ, 
170, 207 

Condon J. J., 1984, ApJ, 287, 461 

Dallacasa D., Bondi M., Della Ceca R., Stanghellini C., 1997, 
Mem Soc astr Ital, 68, 55 

Davis R. J., Unwin S. C., Muxlow T. W. B., 1991, Nature, 321, 
374 

Davis M. M., 1971, AJ, 76, 980 

de Ruiter H. R., Parma P., Fanti C., 1990, AA, 227, 351 
Donnelly R. H., Partridge R. B., Windhorst R. A., 1987, ApJ, 
321, 94 

Doroshkevich A. G., Longair M. S., Zeldovich Y. B., 1970, 
MNRAS, 147, 139 

Dunlop J. S., Peacock J. A., 1990, MNRAS, 247, 19 
Fanaroff B. L., Riley J. M., 1974, MNRAS, 167, 31P 
Fanti C., Fanti R., 1990, in Fanti C., Fanti R., O’Dea C. P., 
Schilizzi R. T., eds, CSS and GPS Radio Sources. CNR - 
Istituto di Radioastronomia Bologna, p. 215 
Fanti C., Vigotti M., Di Paolo C., 1996, in Snellen I. A. G., 
Schilizzi R. T., Rottgering H., Bremer M. N., eds, The 
Second Workshop on GPS and CSS Radio Sources. Leiden 
Observatory, p. 38 

Fomalont E. B., Windhorst R. A., Kristian J. A., 

Kellermann K. I., 1991, AJ, 102, 1258 
Ghisellini G., Padovani P., Celotti A., Maraschi L., 1993, ApJ, 
407, 65 

Giovanini G., Ferreti L., Venturi T., Kom K.-T., Kronberg P. P., 
1993, ApJ, 406, 399 

Giovannini G., Cotton W. D., Feretti L., Lara L., Venturi T., 
Marcaide J. M., 1995, in Cohen M. H., Kellermann K. I., 
eds, Quasars and AGN: High Resolution Imaging. Publ Nat 
Acad Sci, p. 11356 

Gregory P. C., Condon J. J., 1991, ApJ supp, 75, 1011 
Gruppioni G., Zamorani G., de Ruiter H. R., Palma P., 

Mignoli M., Lari C., 1997, MNRAS, 286, 470 
Hales S. E. G., Baldwin J. E., Warner P. J., 1988, MNRAS, 234, 
919 

Hine R. G., Longair M. S., 1979, MNRAS, 188, 111 
Kapahi V. K., Athreya R. M., Subrahmanya C. R., 

McCarthy P. J., van Breugel W., Baker J. C., 

Hunstead R. W., 1996, in Ekers R., Fanti C., Padrielli L., 
eds, Proc IAU Symp 175, Extragalactic Radio Sources. 
Kluwer Academic Pulishers, p. 393 
Kuhr H., Pauliny-Toth I. I. K., Nauber U., 1981, AA Supp Ser, 
45, 367 

Kollgaard R. I., Wardle J. F. C., Roberts D. H., Gabuzda D. C., 
1992, AJ, 104, 1678 

Laing R. A., Peacock J. A., 1980, MNRAS, 190, 903 
Laing R. A., Riley J. M., Longair M. S., 1983, MNRAS, 204, 151 
Laing R. A., Wall J. V., Jenkins C. R., Unger S. W., 1994, in 
Bicknell G. V., Dopita M. A., Quinn P. J., eds, The First 
Stromlo Symposium: The Physics of Active Galaxies. ASP 
Conf Ser 54, p. 201 

Ledlow M. J., Owen F. N., 1996, in Ekers R., Fanti C., 

Padrielli L., eds, Extragalactic Radio Sources. Kluwer 
Academic Publishers, p. 238 
Longair M. S., 1966, MNRAS, 133, 421 
Masson C. R., Wall J. V., 1977, MNRAS, 180, 193 
Metcalfe N., Shanks T., Fong R., Roche N., 1995, MNRAS, 273, 
257 

Moffat A. T., Gubbay J., Robertson D. S., Legg A. J., 1972, in 
Evans D. S., ed, Proc IAU Symp 44: External Galaxies and 
Quasi-Stellar Objects. Reidel: Dordrecht, p. 228 


Morganti R., Killeen N. E. B., Tadhunter C. N., 1993, MNRAS, 
263, 1023 

Morganti R., Oosterloo T. A., Reynolds J. E., Tadhunter C. N., 
Migenes V., 1997, MNRAS, 284, 541 
Morisawa K., Takahara F., 1987, MNRAS, 228, 745 
Murphy D., Browne I. W. A., Per ley R. A., 1993, MNRAS, 264, 
298 

Orr M. J. L., Browne I. W. A., 1982, MNRAS, 200, 1067 
Owen F. N., Ledlow M. J., Keel W. C., 1996, AJ, 111, 53 
Pauliny-Toth I. I. K., Witzel A., Preuss E., Kuhr H., 

Kellermann K. I., Fomalont E. B., Davis M. M., 1978, AJ, 
83, 451 

Pauliny-Toth I. I. K., Steppe H., Witzel A., 1980, AJ, 85, 329 
Peacock J. A., Gull S., 1981, MNRAS, 196, 611 
Peacock J. A., 1987, in Kundt W., ed, Astrophysical Jets and 
their Engines. Reidel: Dordrecht, p. 185 
Perez-Fournon I., Biermann P., 1984, AA, 130, L13 
Press W. H., Teukolsky S. A., Vetterling W. T., Flannery B. P., 
1992, Numerical Recipes in Fortran. Cambridge University 
Press 

Readhead A. C. S., 1995, in Cohen M. H., Kellermann K. I., 
eds, Quasars and AGN: High Resolution Imaging. Publ Nat 
Acad Sci, p. 11447 

Rixon G. T., Wall J. V., Benn C. R., 1991, MNRAS, 251, 243 
Robertson J. G., 1980, MNRAS, 190, 143 

Roche N., Shanks T., Metcalfe N., Fong R., 1996, MNRAS, 280, 
397 

Rowan-Robinson M., Benn C. R., Lawrence A., 

McMahon R. G., Broadhurst T. J., 1993, MNRAS, 263, 123 
Sadler E., Jenkins C. R., Kotanyi C. G., 1989, MNRAS, 240, 591 
Saunders W., Rowan-Robinson M., Lawrence A., Efstathiou G., 
Kaiser N., Ellis R. S., Frenk C. S., 1990, MNRAS, 242, 318 
Scheuer P. A. G., Readhead A. C. S., 1979, Nature, 277, 182 
Scheuer P. A. G., 1987, in Zensus J. A., Pearson T. J., eds, 
Superluminal Radio Sources. Cambridge University Press, 
p. 104 

Schmidt M., 1976, ApJ, 209, L55 

Shaver P. A., Wall J. V., Kellermann K. I., Jackson C. A., 
Hawkins M. R. S., 1996, Nature, 384, 439 
Shaver P. A., Hook I., Jackson C. A., Wall J. V., 

Kellermann K. I., 1998, in Carilli C., Radford S., Menten K., 
eds, Highly Redshifted Radio Lines. ASP Conf Ser, in press 
Singal A. K., 1996, MNRAS, 278, 1069 
Snellen I. A. G., 1997, PhD thesis, University of Leiden 
Taylor G. B., Vermeulen R. C., Pearson T. J., 1995, in 

Cohen M. H., Kellermann K. I., eds, Quasars and AGN: 
High Resolution Imaging. Publ Nat Acad Sci, p. 11381 
Urry C. M., Padovani P., 1995, PASP, 107, 803 
Vermeulen R. C., 1995, in Cohen M. H., Kellermann K. I., eds, 
Quasars and AGN: High Resolution Imaging. Publ Nat 
Acad Sci, p. 0 

Wall J. V., Jackson C. A., 1997, MNRAS, 290, 17P 
Wall J. V., Peacock J. A., 1985, MNRAS, 216, 273 
Wall J. V., Pearson T. J., Longair M. S., 1980, MNRAS, 193, 
683 

Wall J. V., 1980, Phil Trans R Soc, A296, 367 
Wall J. V., 1994, Aust J Phys, 47, 625 
Wall J. V., 1998, in Jackson N., Bremer M. N., 

Perez-Fournon I., eds, Observational Cosmology with the 
New Radio Surveys. Kluwer Academic Publishers, p. 129 
Windhorst R. A., Miley G. K., Owen F. N., Kron R. G., 

Koo D. C., 1985, ApJ, 289, 494 
Wrobel J. M., Krause S. W., 1990, ApJ, 363, 11 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science DT^X style file. 


© 1998 RAS, MNRAS 000, 0L1 


